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Abstract 

The interaction of the atmosphere with the solar radiation 
reflected from the Earth's surface heavily affects the remotely 
sensed images acquired by passive optical sensors. In order 
to obtain unbiased estimations of the surface reflectance, the 
above image degradation and radiometric distortion of the 
acquired data must be corrected, an operation which is called 
"atmospheric correction" often. Knowing the local 
parameters of the Earth atmosphere is therefore necessary for 
performing any atmospheric correction procedure, moreover 
this knowledge also is important for many environmental 
studies and remote sensing applications. In this paper a 
direct method of estimating the atmospheric optical depth 
(direct transmittance) using images acquired by the sensor 
CHRIS is discussed. The estimation of the optical depth is 
possible since this instrument is able to acquire up to five 
overlapping images of the same ground location observed 
under different sightlines. In such a way, this method allows 
us to estimate the atmospheric optical depth for twelve 
spectral channels of the CHRIS spectrometer. Such estimates 
can be used either for an optimal tuning of the atmospheric 
correction algorithm or to assess the state of the atmosphere. 

Keywords 

Component; Optical Depth; Atmospheric Correction; Atmospheric 
Optics; Remote Sensing; CHRIS-PROBA 

Introduction 

For a long time the remotely sensed data acquired by 
passive optical sensors have been used for answering 
many questions in a wide range of disciplines that are 
concerned with land-surface processes (Asrar 1989). 
Such instruments are known to obtain Earth's surface 
images in several narrow bands in visible and infrared 
regions of the electromagnetic spectrum. This 
information allows us to state and solve the problem of 
the Earth's surface identification. The accuracy of such 
deductions is determined by many factors depending 
on both sensor characteristics and conditions of 
observation. One of the most important factors is the 



influence of the atmosphere on the radiation collected 
by sensor. The interaction of the solar radiation 
reflected from the ground with the atmosphere hinders 
the quantitative use of the remotely sensed data, 
demanding for the atmospheric correction of the 
collected data. In its turn, this correction requires the 
knowledge of various atmospheric and signal 
parameters, which concur to determine the signal 
distortion caused by the atmosphere. The most 
straightforward solution to this problem is obtained 
when these parameters can be directly deduced from 
the remotely sensed data alone. In this case the 
interfering effect on the atmospheric parameter 
estimation originated by reflection / scattering 
phenomena taking place near the ground may be very 
strong. Nevertheless, in certain specific cases some 
atmospheric characteristics can be measured 
(estimated) from remotely sensed data. A successful 
example of such estimation is the atmospheric 
correction of the ocean color, where the very low 
reflectance of the water in the red region of the visible 
spectral range can be used for the atmospheric 
correction (Gordon 1978). A number of papers and 
research works have been devoted to the estimation of 
atmospheric optical parameters and the atmospheric 
correction of remotely sensed images (Liang et al. 2001, 
Kaufman 1984, Lyapustin 2001, Thome et al. 1998, 
Kaufman et al. 1997, Davis et al. 2010, Karpouzli and 
Malthus 2003, Wen et al. 1999, Barducci et al. 2008, 
Semenov et al. 2011). 

The estimation of the atmospheric characteristics 
relying on inverse modeling of remotely sensed data 
only is the preferable practice because it doesn't 
require any additional equipment, and permits to 
recover correct information from past acquisitions. 
However, it should be emphasized that these estimates 
can be obtained only for images having specific 
properties in terms of observed materials and sensor's 
characteristics, and that the accuracy of such estimates 
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cannot be reliably controlled always. Let us note that 
most remote sensing instruments are designed for 
Earth observation, not for the measurement of 
atmospheric features. Nevertheless, such estimates are 
a possible source of additional and sometime unique 
information about the local state of the atmosphere. 

In this paper we discuss a direct method for assessing 
the atmospheric optical depth, i.e. the direct 
transmittance, which adopts the remotely sensed data 
only. The paper is organized as follows. Section II 
introduces the basic method for estimating the 
atmospheric optical depth stemming from remotely 
sensed images acquired under two angles of 
observation, and shows its limitation. The extension of 
this method to images acquired from five view points, 
as in a case of CHRIS/PROBA is also addressed. 
Section III shows the results of atmospheric optical 
depth estimation obtained from natural images in 12 
spectral channels of the imaging spectrometer CHRIS- 
PROBA. The last section draws some conclusions. 

Estimation of Atmospheric Optical Depth 

Our method of estimating the atmospheric optical 
depth t , or the direct transmittance T=e~ T , is based 
on the comparison between the two values of the 
radiance variance averaged over all the pixels of the 
same area in the two images received at two different 
observation angles (Semenov 2002, Moshkov et al. 
2009). 

Estimation of Atmospheric Optical Depth Using 
Images Acquired from Two Observation Angles 

Suppose we have two images of the same ground 
location acquired under the observation angles 3 l and 

3 2 in the same spectral channel of a passive optical 
sensor. Without loss of generality we can assume 
3 l = and 3 2 = 3 , so one image corresponds to the 
nadir viewing geometry. 

The radiance of an arbitrarily resolved pixel of the 
image gathered over the nadir direction can be written 
as 

b^=b^+B^e- T i =1,2... N 

where by is the atmospheric path radiance not 
reflected from the surface, is the radiance of the ; th 
resolved element on the surface observed in nadir, and 
N is the number of pixels in the image. Here, we 
assume the image have been preliminary corrected to 
remove the scattering contribution to radiance (the so- 



called "adjacency effect", Semenov et al. 2011). 
Moreover, we assume that the image size (i.e. the 
sensor FOV) is small enough that possible variations of 
the optical depth along the sightline corresponding to 
different pixels can be disregarded. 

When observing the same resolved ground location 
from the slant sightline 3 we obtain the radiance 
stated by the following equation: 

bW=b<p + BWe->" i = \,2,-N 

where b ( f and are the atmospheric path 

radiance, and the radiance emitted by the resolved 
element observed at the angle 3 , the air-mass 
parameter being represented by ju = 1/cos 3 . 

The variance of the radiance image seen in the nadir 
direction can be determined as 

where m^= — X^,- ' is the average for the nadir 
N i=i 

image. Similarly, the radiance variance for the image at 
the angle 3 can be calculated as: 

s(*) = ±2tf*)- m w? 

V ' N i=i L J 

where m^= — Y.b\ 9 ^ is the image average 

N i=i 

corresponding to the observation angle 3 . The 
estimated atmospheric optical depth f can then be 
defined as shown in the next equation: 

2k(3) {s(0)j 

And k(3)=ju-l. 

Let's now consider the errors of the estimate (5). 
Evidently, the radiance of an arbitrarily resolved 
element of the image at the angle 3 can be written as: 

bf = b^ + [B^ + p^]e-^, i=l,2,...N 

where fyf' is the difference of emitted radiance for the 
z-th pixel between the two viewing directions 3 and 0. 
The image average of the radiances coming from the 
two datasets (view-lines) yields: 

m^=b^+M^e- T 
m^=b[ 3 KM^e-^ + AM^e-^ 

where M^= — 2Z?(°) is the mean emitted radiance 
N i= i 

for the image observed in the nadir direction 
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(neglecting the atmospheric attenuation), and 
AM w = — Y. 0r^ is the radiance deviation averaged 



N i= i 



over the image at the angle 3 . From (7) it is possible to 
compute the deviation Abj^ of the radiance emitted by 
a generic pixel from its image average, as shown in the 
equation below: 



Ab\ 
Abi 



b /o)_ m (o) : 



m 



AB 



+ Ap\ 



where AB ( ( ) = -M® and A#=#-AM^. In 
view of (8), it is easy finding an expression for the 
radiance variance S(&) : 



It 



N i= \ 



V ; N U 

Then the estimate of the atmospheric optical depth r 
in (5) becomes: 



**) = - 



-In 



S(S) 



-T + n r 



where 



2k{9) 



In 



2k{3) {S{0) 

2Zab( )a^) + z[a^)] 2 



1 + 



i=l 



is the estimation error, and r is the true value of the 
atmospheric optical depth. Equation (11) just implies 
that in the case of a Lambertian field of reflectance 
( A/3^ = ) the estimate (10) equates the true value of 
the atmospheric optical depth, while any deviation of 
the reflectance of the observed scene from the behavior 
of a Lambertian surface gives rise to possible errors 
affecting i . It should also be noted that equations 
from (5) to (11) have been obtained assuming the two 
images involved in calculations are exactly matched 
(co-registered), otherwise an additional error due to 
spatial mismatching arises. 

The following conclusions can be derived from the 
analyses of the eq. (11). 

1) For surfaces with fixed scattering properties 
the lesser is k(&) (the lesser is an angle of 
observation <9 ) the greater is the estimation 
error n T . 

2) The error n T only depends on the centered 
value of the involved reflectance fields. 



Therefore, if we assume AS,( 0) and ApW to be 
independent random variables, the sum 

in (11) is zero, and the greater is 

the amplitude of AB,( o) with respect to Afi\ 9 \ 
the lesser is the error. For a given Afi\ & \ an 
augment of ABf^ corresponds to an increase of 
the range of spatial variability along the image, 
i.e. increasing the variance of the image. 

3) Since the angular scattering properties of the 
natural surfaces constituting the available 
images are unknown, we suppose that the 
estimation of the atmospheric optical depth r 
deduced from two different sight lines (images) 
only is prone to large errors. 

Estimation of Atmospheric Optical Depth Using 
CHRIS Images (Multiple Observation Angle) 

A new generation of optical sensors has been 
developed in the last years. Such sensors are 
characterized by a large number of spectral channels, 
finer radiometric accuracy, lesser volume and mass 
budgets, and in some cases the ability to capture 
Earth's images of the same site for various angles of 
observation. A good example of such sensors is the 
Compact High Resolution Imaging Spectrometer 
(CHRIS) developed by the European Space Agency 
(ESA), and mounted onboard the compact satellite 
PROBA. CHRIS-PROBA operates over the visible and 
near infrared spectral region, from 400 nm up to 1050 
nm. This spectral range is covered with 63 spectral 
bands at a Ground Sampling Distance (GSD) of 36 m, 
or with 18 bands at a spatial resolution of 18 m. The 
actual CHRIS spatial resolution is somewhat variable 
as its altitude varies from 553 km to 588 km. Using 
PROBA's angle steering capabilities in along and 
across track directions, sets of 5 images can be acquired 
for the selected site along the same orbit, at the 
following Fly-by Zenith Angles (FZA, along track) of 
+55.0°, +36.0°, 0.0°, -36.0°, -55.0°. Additional information 
regarding the CHRIS-PROBA characteristics and its 
acquisition modes can be found in Barducci et al. 2009. 
In order to adopt CHRIS-PROBA data for autonomous 
estimation of the atmospheric optical depth, we have 
further developed the optical depth estimation 
procedure. The resulting modeling is free from the 
limitations characteristic of the basic procedure 
described above that only uses images from two 
observation angles. 
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Let us suppose to have five matched images of the 
same area obtained in the same spectral channel for the 
five observation geometries typical of the CHRIS- 
PROBA sensor, indicated by angles 3 =0.0, 3 t =36.0, 
3 2 =55.0, #3 =-36.0, and 3 4 =-55.0 . Hence, we can 
compute four independent estimates of the 
atmospheric optical depth for the four non-null 
observation angles: 

The surfaces of targets framed in the available images 
will not have an exact Lambertian behavior, thus the 
various estimates resulting from (12) will be 

mutually different, and each one of them will differ 
from the true value r . As long as the statistical 
properties of the measurement "noise" are unknown, it 
is appropriate to accept as optimal r estimation the 
value that minimizes the dispersion of the estimates 
(the standard deviation of the estimates). As known, 
the arithmetic mean is the estimator that obeys this 
property. However, considering the dependence 
reported in (11) of the estimation error on the 
observation angle, it seems desirable to adopt a 
weighted average of the estimates obtained in various 
viewing geometries. To this purpose, we consider the 
general weighted dispersion formula: 

N (- - \2 
I.C0j{Tj-T m ) 

7=1 

that is minimized by the following estimator: 

7Y 

2>7 f 7 

T 

"m N 

7=1 

In view of (11) we choose the weights «. = k 2 ($j), 
which have the property k(3 l ) = k(& 3 ) and k{3 2 ) = k{3 A ). 
In summing up, we can write: 

# _ k 2 {3 l ){r l +T i )+k 2 {3 2 \r 2 +i A ) 
2[fe 2 (^)+fe 2 (f 2 )] 

Equation (15) holds true for an atmospheric optical 
depth retrieved from CHRIS-PROBA data, and is easily 
extended to the general case of N available viewing 
directions, i.e. (14). It is worth noting that the optical 
depth estimation involving data from multiple Point of 
Views (POVs) allows us to increase the accuracy of the 
final estimation, and even to control its quality. When 
observing an ideal Lambertian surface, any estimation 



matches the actual value of r , and in the 

opposite circumstance i(&j)=t r(<9/) Y/'#Z . Therefore, 

the dispersion of the N available estimations is a 
measure of the quality of the average estimate in (14- 
15). We can define the following discrepancy function 
S which measures the quality of the inferred optical 
depth: 

7=1 

Let us note that a vanishing s is a necessary condition 
of an exact estimate i [3 ■ ) = t Vy* , but it is not sufficient. 

This means that t (i9 ■ ) = r V/ implies s = , but the 

converse is not necessarily true. For the specific case of 
the CHRIS-PROBA images, for which 3 X = -3 3 and 
3 2 = ~3 4 , the estimation quality can also be 
characterized by the next expression that defines the 
discrepancy of means e p : 

7=1 

t mJ ={f j+ t j+2 )/2 j =1,2 

Unfortunately, as in a case of the ordinary discrepancy, 
zero discrepancy of means is a necessary but not at all 
sufficient condition for the equality between the 
estimate i m and the true value r . Let us note that e p 

can be zero even for a non-lambertian surface, when 
t ml =r, 1 = 1,2- 

Since Lambert reflection is an idealization of a complex 
phenomenon, only rarely this behavior will account for 
all the targets of a whole scene. At the same time, most 
images contain a great number of different objects, 
both natural and artificial, with various scattering 
(reflecting) properties. Repeating the atmospheric 
optical depth estimation for a number of single image 
fragments significantly increases the probability to find 
an image fragment with quasi-Lambertian reflectance. 
The estimation performed for this fragment should be 
characterized by null or vanishing discrepancy, and we 
could select the minimal discrepancy estimation from 
the mess of estimations obtained in this way. 

Optical Depth Estimation for a CHRIS 
Dataset 

The optical depth estimation procedure previously 
outlined has been applied to a CHRIS-PROBA dataset 
acquired over the San Rossore test site (Pisa, Italy), on 
the 16th of June 2003. For this acquisition, the CHRIS 
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spectrometer was operated in Mode 4, collecting data 
from 18 spectral bands at a GSD of 18 m. Unfortunately, 
the last six spectral channels of this and other CHRIS- 
PROBA acquisitions experienced significant calibration 
troubles (Barducci et al. 2009), a circumstance that 
hindered the option to retrieve significant optical 
depth estimations at these wavelengths. Validation 
data collected in-field gave the following results: 
temperature +28.50 C, relative humidity 65%, pressure 
1011 mbar, integrated solar irradiance 810 w nr 2 , and 
diffuse solar irradiance 170 w m 2 (both from 0.285 pm 
- 2.800 pm). Fig. 1 shows a true color composite of a 
mosaic of the five angular views provided by the 
CHRIS-PROBA imaging spectrometer. 




Fig. 1 Example of CHRIS-PROBA acquisition over the San 
Rossore (Pisa, Italy) test site. The central image (top of picture) 
has been gathered with zero fza, the two images on the 
central ragion of the picture have fza=±36° and the two 
bottom images correspond to fza=±55°. images on the left 
column have positive fza; e.g. the image in the bottom left 
corner is the san rossore view obtained with fza=+55°. 

We have estimated the atmospheric optical depth for 
12 CHRIS spectral channels (Moshkov et al. 2009), 



without applying compensation for the adjacency 
effect. The optical depth estimation procedure 
consisted of two phases: co-registration of the five 
images constituting the dataset, and the procedure of 
the optical depth estimation discussed above. 

In the co-registration phase we have adopted the affine 
transformation shown in the following equations: 

x = ax + by + c 
y' = dx + ey + f 

Here x and y are the coordinates of a pixel in the first 
image, and x' and y' are those in the second one. 
Coordinates are, as usual, column and row indexes 
( x, y ), as measured from the top left corner of each 

image. In order to find the coefficients in (18), some 
Ground Control Points (GCPs) have been visually 
chosen in both images, and their coordinates employed 
to calculate the optimal transformation connecting the 
two maps (we used eight points spread over the entire 
image). This yields an over-determined system of 16 
equations with only six unknowns whose solution can 
be inferred by the least-squares method. With this 
procedure we have found the transformation 
coefficients (a,b,...,f ) from the nadir image to each of 
the four slant path images completing the dataset. The 
least squares minimization has been iterated in order 
to increase the precision of the co-registration phase. 
The accuracy of the determination of the 
transformation coefficients has been assessed as the 
r.m.s. distance between measured and predicted GCP 
coordinates. 

The structure of our algorithm for r estimation 
reminds the moving window procedure often adopted 
for processing digital images. Our algorithm can be 
applied to monochromatic multi-angular datasets, 
meaning that multiple band datasets are processed 
sequentially iterating the algorithm for each spectral 
channel. The sequence of operations constituting the 
application of the algorithm to a single spectral band of 
a multi-angular dataset are depicted in the following. 
A square window, with size N p. pixel, is adopted for 

scanning the available nadir image, with unit step 
along both axes. Window sizes in the range from 84 
pixels up to 116 pixels have been attempted. The 
scanning procedure of a monochromatic image 
incorporates the following operations. 

1. Extracting from any image of the multi-angular 
set the sub images corresponding to the current 
position of the moving window. This operation 
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utilizes the affine transformation in (18) with 
the free parameters previously determined in 
the co-registration phase. 

2. Computing the variance in each multi-angular 
view of pixels conjugated to positions inside the 
moving window. 

3. Estimating the optical depth f[&j) for any 

available observation angle (12), the optimal 
optical depth f m (15), the corresponding 
discrepancy s (16), and discrepancy of means 
s p (17). 

4. Steps from 1 to 3 are iterated after changing the 
position of the moving window. To reduce the 
computation burden we employed an empirical 
variance threshold a in order to prevent the 
optical depth calculation over image fragments 
of lesser dispersion (which would produce 
lower quality optical depth estimates). At the 
completion of the scanning process associated 
to a given window size, the optical depth 
estimation i ml corresponding to the minimal 
discrepancy s { = min{^} is retained together 
with the discrepancy s t itself, and the 
discrepancy of means e t . 

5. Once the image scanning is completed, the 
entire procedure, steps from 1 to 4, is repeated 
using a larger moving window size. Window 
size has increased with a step of 4 pixels. 

This originates nine optimal estimates f mI , nine 

minimal discrepancies £ t , and nine discrepancies of 

means £ t , each one connected to the nine admitted 

sizes of the moving window. Among them we select as 
the best estimate i mm of the optical depth that optimal 

estimate T ml corresponding to the minimal 

discrepancy of means £ pl , as shown by the next 

equation: 

tmm= / = 1 9 ) 

The optical depth estimates r mm and the 
corresponding values of the direct transmittance 

T = e~ f ' obtained from the CHRIS-PROBA images in 

12 spectral channels are detailed in Table I. The 
corresponding atmospheric direct transmittance 
simulated for the same twelve spectral channels by 
means of the MODTRAN 5.2 (Anderson et al. 1995) 



code are plotted in Fig. 2 together with the 
transmittance estimates. Simulations have been 
executed for the Mid-latitude summer model, utilizing 
four different aerosol conditions: 1) rural, with 
visibility of 23 km, 2) marine aerosol with 23 km of 
visibility, 3) rural aerosol and visibility of 5 km, and 4) 
urban aerosol with 5 km of visibility. Optical depth 
estimation in the last six spectral channels gave 
evidently wrong results, presumably due to the 
ascertained degradation of the sensor's radiometric 
response affecting these bands, see Barducci et al. (2009) 
and Cutter (2004). Let us note that the peak and the 
average Signal-to-Noise Ratio (SNR) typical of CHRIS- 
PROBA datasets decrease with increasing the 
wavelength (Barducci et al. 2005), giving rise to an 
additional possible reason for explaining the poor 
estimations obtained in the near infrared bands. 



Table I Optical depth estimation from CHRIS-PROBA images 





Km ( nm ) 


T 

* mm 


f 


T (Modtran, 15 
km visibility) 


1 


489.2 


0.736 


0.479 


0.498 


2 


551.2 


0.594 


0.552 


0.549 


3 


631.2 


0.516 


0.597 


0.602 


4 


671.7 


0.473 


0.623 


0.653 


5 


679.9 


0.452 


0.636 


0.659 


6 


685.6 


0.459 


0.632 


0.628 


7 


691.3 


0.468 


0.626 


0.607 


8 


697.2 


0.420 


0.657 


0.642 


9 


703.2 


0.396 


0.673 


0.634 


10 


709.3 


0.362 


0.696 


0.661 


11 


715.5 


0.397 


0.672 


0.621 


12 


734.9 


0.302 


0.739 


0.649 
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Fig. 2 Example of atmospheric optical depth estimation from 
chris-proba data, and its comparison with several modtran 
5.2 simulations calculated for various atmospheric models 
and several visibility values. 
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As can be seen the estimates compare fairly with the 
simulations, reproducing a similar amplitude and 
spectral distribution. The transmittance estimates 
deduced from CHRIS-PROBA data nearly approximate 
the Modtran 5.2 simulations with 15 km of visibility, as 
can be seen in Fig. 2 and Table I. 

Conclusion 

The Earth's atmosphere heavily affects the remote 
sensing images acquired by spaceborne passive optical 
sensors due to radiation-matter interaction phenomena 
like radiation absorption and scattering. One of such 
phenomena is the direct attenuation of radiation 
reflected from the Earth's surface when it's passing 
through the atmosphere. The attenuation is 
characterized by the optical depth, which in turn 
determines the direct transmittance of the atmosphere. 

In this paper we have suggested a new method for 
estimating the atmospheric optical depth stemming 
from multiangle remotely sensed images, not using 
any in-field measurement or external information 
about the atmosphere and the Earth's surface. We have 
reviewed possible errors affecting the optical depth 
estimations obtained with this method. The realization 
of this method permitted us to estimate the optical 
depth (direct transmittance) for twelve spectral 
channels of a data volume gathered by the imaging 
spectrometer CHRIS-PROBA. Estimations as far 
obtained compare fairly with optical depth simulations 
computed by the MODTRAN 5.2 code, especially for 
the first CHRIS-PROBA spectral channels. 

We remark that optical depth values directly extracted 
from remotely sensed hyperspectral and multiangular 
images are a supplementary source of information 
concerning the state of the atmosphere. Potentially, the 
availability of such estimations fills the gap between 
the data acquisition and their utilization, permitting us 
to mitigate the complexity that hinders the generation 
of high-level remote sensing products. The assessment 
of high-level remote sensing products is often 
unmanageable due to the absence of validation and 
calibration in-field measurements, a shortage that is 
proficiently overcome by the discussed estimation 
algorithm. 
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